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■ Abstract: We study in detail the ability of the nominal configuration of the IceCube neu- 

trino telescope (with 80 strings) to probe the parameter space of the Constrained MSSM 
^ ■ (CMSSM) favoured by current collider and cosmological data. Adopting conservative as- 

^\ . sumptions about the galactic halo model and the expected experiment performance, we 

find that IceCube has a probability between 2% and 12% of achieving a 5a detection of 
dark matter annihilation in the Sun, depending on the choice of priors for the scalar and 
^ ' gaugino masses and on the astrophysical assumptions. We identify the most important 

' annihilation channels in the CMSSM parameter space favoured by current constraints, and 

we demonstrate that assuming that the signal is dominated by a single annihilation channel 
can lead to large systematic errors in the inferred WIMP annihilation cross section. We 
demonstrate that ^ 66% of the CMSSM parameter space violates the equilibrium con- 
dition between capture and annihilation in the center of the Sun. By cross-correlating 
our predictions with direct detection methods, we conclude that if IceCube does detect a 
neutrino flux from the Sun at high significance while direct detection experiments do not 
find a signal above a spin-independent cross section > 7 x lO"*^ pb, the CMSSM will 
be strongly disfavoured, given standard astrophysical assumptions for the WIMP distribu- 
tion. This result is robust with respect to a change of priors. We argue that the proposed 
low-energy DeepCore extension of IceCube will be an ideal instrument to focus on relevant 
CMSSM areas of parameter space. 
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1. Introduction 



Supersymmetry (SUSY) is the most promising candidate for physics beyond the Standard 
Model (SM). It provides solutions to two of the most challenging questions facing particle 
physics today Q. Namely, the instability of the Higgs mass against radiative corrections 
(known as the "fine-tuning problem") and the nature of dark matter. Firstly, the fine- 
tuning problem is addressed via the cancellation of quadratic divergences in the radiative 
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corrections to the Higgs mass. Secondly, the theory predicts the existence of Weakly In- 
teracting Massive Particles (WIMPs). Assuming R-parity conservation, the lightest super- 
symmetric particle is stable and can have survived as a relic from the Big Bang, meeting 
the necessary requirements to account for the the dark matter content of the Universe 
suggested by observations |^]. 

Under the above mentioned assumptions, our galaxy could be embedded in a dark 
halo formed of WIMPs. During the last decades a wide range of experiments have been 
developed with the goal of detecting such WIMPs. The overall strategy has followed a 
double-track approach: on the one hand, "direct detection experiments" aim at measuring 
the energy deposited by the interaction of a WIMP from the halo with a nucleus in the 
target material |^]. On the other hand "indirect detection experiments" focus on the 
detection of the annihilation products of WIMPs, such as cosmic rays, neutrinos or 7 rays, 
produced by dark matter annihilation from a variety of sources, including the galactic 
center satellite dwarf galaxies |^], and diffuse emission from unresolved sources. 

A promising search technique is to look for an excess of neutrinos from WIMP an- 
nihilations in the center of heavy bodies, like the Earth or the Sun. The idea behind 
this possibility is that WIMPs from the halo can be gravitationally trapped in orbits in 
the Solar System and, by scattering off the matter of heavy bodies, fall below the escape 
velocity and accumulate in their center Q. These trapped WIMPs can then annihilate 
with each other to produce quarks, leptons, gauge or Higgs bosons. The decays of some 
of these particles will produce neutrinos, which will escape the body, providing an excel- 
lent experimental signature: they will be directional and will have energies from tens to 
several hundreds GeV depending mainly on the WIMP mass. In the absence of a signal, 
the existing neutrino telescopes, AMANDA, IceCube, ANTARES and Baikal, as well as 
the SuperKamiokande detector have already set limits on the neutrino flux from neutralino 
annihilations in the Sun or Earth [0, |^ |^, [l^, 

As mentioned above, the supersymmetric framework offers a natural and well moti- 
vated candidate for a WIMP, namely the neutralino. The neutralino is the mass eigenstate 
of a mixture of bino, wino (the superpartners of the B, gauge bosons) and higgsinos 
(the superpartners of the H^, H2 Higgs bosons). Whenever the lightest among the four 
neutralinos is the lightest SUSY particle, it is expected to be stable due to R-parity (in the 
following, "neutralino" always denotes the lightest of the 4 neutralino eigenstates). There- 
fore the neutralino provides a stable, neutral, massive (typical masses range from tens of 
GeV up to several TeV) particle that can account for the cold dark matter required to fit 
cosmological observations. 

Despite its potentiality, low-energy SUSY has a fundamental question yet to be an- 
swered; namely, how it is broken. Along the past years, several SUSY-breaking mechanisms 
have been proposed in the literature. Perhaps the most appealing one regards to gravity 
as the mean of communication among a hidden sector where SUSY breaking occurs and 



the visible sector [12[. This scenario is commonly called Supergravity, and a simple realiza- 



tion of it is the so called Constrained Minimal Supersymmetric Standard Model (CMSSM) 



- 2 - 



fl^ , It is defined in terms of five free parameters: common scalar {itiq), gaugino 

(mi/2) and tri-linear (Aq) mass parameters (all specified at the GUT scale) plus the ra- 
tio of Higgs vacuum expectation values tan/3 and sign(/i), where /i is the Higgs/higgsino 
mass parameter whose square is computed from the conditions of radiative electroweak 
symmetry breaking. 

There is a vast literature devoted to the study of the phenomenology of the CMSSM, 
based on current collider results and data from dark matter searches. The usual proce- 
dure to explore the model's parameter space has been to evaluate the likelihood of the 



observables on a a fixed grid, often encompassing only 2 or 3 dimensions at the time |15|. 
This method has several limitations. Firstly, the number of points required scales as , 
where is the number of the model's parameters and k the number of points for each 
of them. Therefore this technique becomes highly inefficient for exploring with sufficient 
resolution parameter spaces of even modest dimensionality. Secondly, some features of the 
parameter space can easily be missed by not setting a fine enough resolution. Thirdly, extra 
sources of uncertainties {e.g., those due to the lack of precise knowledge of SM parameter 
values) and relevant external information {e.g., about the parameter range) are difficult to 
accommodate. 

More recently a new approach based on Bayesian methods (for a recent review of 
Bayesian methods, see, e.g. [^) has been successfully applied to the exploration of su- 
persymmetry phenomenology 0, 111, pi, M M, l2|, |1, pi l2|, g pi l2|, IsO, |T|. 



Bayesian methods coupled with Markov Chain Monte Carlo (MCMC) technology are su- 
perior in many respects to traditional grid scans of the parameter space. Firstly, they 
are much more efficient, in that the computational effort required to explore a parameter 
space of dimension scales roughly proportionally with N for standard MCMC meth- 



ods. The more recent "nested sampling" algorithm 33 1 offers an efficiency which is 
roughly constant with N, thus opening up the possibility of exploring parameter spaces 
of high dimensionality such as the general MSSM. Secondly, the Bayesian approach allows 
one to incorporate easily into the final inference all relevant sources of uncertainty. For 
a given SUSY model one can include relevant SM (nuisance) parameters and their asso- 
ciated experimental errors, with the uncertainties automatically propagated to give the 
final uncertainty on the SUSY parameters of interest. Thirdly, another key advantage is 
the possibility to marginalize (i.e., integrate over) additional ("hidden") dimensions in the 
parameter space of interest with essentially no computational effort. Finally, a key feature 
is the ability to obtain probability distributions for any function of the input parame- 
ters, which makes it possible to generate quantitative statistical predictions within a given 
SUSY model for any observable quantity of interest. This most useful feature is exploited 
in this paper to produce statistical predictions for the capability of IceCube to constrain 
the CMSSM. On the other hand, the priors issue, considered as a caveat of the Bayesian 



approach, has now been thoroughly addressed in this context |2^, 34, 36]. 



Previous works have considered the detection prospects for neutrino telescopes (e.g., 
1 37 , |3^, ^ ) , but without adopting a statistical framework allowing to quantify the prob- 



ability of discovery within a given theoretical model. More recently, Refs. p^, E0| have 
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employed a similar technique as the one proposed here, but presented their predictions in 
terms of muon flux at the detector location. There is a vast literature concerning indirect 
detection signatures of neutrino from the Sun, including models beyond the CMSSM, see 
e.g. |^T| (we refer to and references therein for a more exhaustive list of works on the 
subject). In this study we improve on previous works by focusing our analysis on the 
number of muon events in the detector, and by showing that this procedure is immune to 
potentially large systematic errors that can affect the comparison between theory and data 
when the dominant annihilation channel is unknown. 

The paper is organized as follows. In section ^ we introduce our scanning methodology 
and our statistical approach, including the data used to constrain the CMSSM parameter 
space. We then review the basic elements about indirect WIMP searches using neutrino 
telescopes (section ^) . We present our findings about indirect detection prospects for Ice- 
Cube and analyze in detail the contribution to the expected signal from several annihilation 
channels (section ^). We compare the indirect detection prospects with direct detection 
methods in section |5| and we discuss the dependency of our results on priors assumptions 
in section |^. We present our conclusions in section |7[ 



2. Methodology 



2.1 Statistical framework and scanning method 

We denote the parameter set of the CMSSM introduced above, (mo, 5 and tan /?), by 
9 (we fix sgn(;u) to be positive, motivated by arguments of consistency with the measured 
anomalous magnetic moment of the muon), while ijj denotes the relevant SM quantities 
that enter in the calculation of the observable quantities (so-called "nuisance parameters" ) , 
namely 

V^ = {Mi,m,(m,)^,a,(Mz)^,aem(Mz)^}, (2.1) 

where Mt is the pole top quark mass, mb[mi,)^^^ is the bottom quark mass at mi,, while 
OemiMz)'^^ and as{Mz)^^ are the electromagnetic and the strong coupling constants 
at the Z pole mass Mz, the last three evaluated in the MS scheme. We denote the full 
8-dimensional set of parameters by 

G = (0,^). (2.2) 
The cornerstone of Bayesian inference is Bayes' Theorem, which reads 

,(e|.).«5). 



The quantity p{Q\d) on the l.h.s. of eq. ( |2.3| ) is called the posterior, while on the r.h.s., the 
quantity p{d\Q) is the likelihood (when taken as a function Q for fixed data, d). In this 



paper, we use the same likelihood function as given in Sec. 3.1 of ref. |19], with updated 
values for the experimental constraints as applicable (see Table ||). The quantity p{Q) is 
the prior which encodes our state of knowledge about the values of the parameters Q before 
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we see the data. The state of knowledge is then updated to the posterior via the hkehhood. 
Finally, the quantity in the denominator is called evidence or model likelihood. If one is 
interested in constraining the model's parameters, the evidence is merely a normalization 
constant, independent of 0, and can therefore be dropped. 



In order to explore efficiently the posterior of Eq. ( ^.31) , we adopt the MultiNest |33| 
algorithm as implemented in the SuperBayeS code [^]. MultiNest provides an extremely 
efficient sampler even for likelihood functions defined over a parameter space of large di- 
mensionality with a very complex structure. This aspect is very important for the CMSSM, 
as previous MCMC scans have revealed that the 8-dimensional likelihood surface is very 
fragmented and that it presents many finely tuned regions that are difficult to explore 
with conventional MCMC (and almost impossible to find with conventional grid scans). 
By explicitly including the SM parameters and then integrating over them at the end, our 
scans automatically account for the full uncertainty associated with the current experi- 
mental errors on those quantities. In particular it has been shown that it is important to 
include the impact of the uncertainty in the top mass to compute accurately the CDM relic 
abundance (see Figure 3 in |p2||). 



2.2 Choice of priors 



In order to perform the scan over the CMSSM and SM parameters, we need to specify a 
prior for the l.h.s of Bayes' theorem, Eq. (|2.3| ). The role of the prior is to define a statistical 
measure on the parameter space. In principle, when the likelihood is strongly constraining 
(i.e., for accurate data) the posterior is dominated by the likelihood and the prior choice 
is irrelevant, as the information in the likelihood completely overrides the information in 
the prior (for an illustration, see Fig. 2 in fl^). However, it has been shown that this 
is presently not the case for the CMSSM, i.e. different (plausible) choices of priors do 
lead to different posteriors and hence different inferences on the parameter space (for a 
detailed analysis, see |Q). The resolution of this issue will come from future data, and in 
particular LHC measurements of the SUSY mass spectrum, which will conclusively resolve 
the ambiguity brought about by the choice of priors |42]. Whenever several prior choices 
appear plausible based on physical considerations, it is important to verify to which extent 
one's statistical conclusions depend upon the choice of priors. In this work, we present our 
main results for the choice of so-called "log priors", namely priors that are flat in logmx/2 
and log rriQ , and flat in tan f3 and Aq . The prior ranges are chosen as to cover the region 
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where the hkehhood is significantly non-zero, and span the fohowing intervals^ : 

(corresponding to 50 MeV < mo,mi/2 < 4 TeV) (2.5) 
- 7 TeV < ^0 < 7 TeV (2.6) 
2<tan/3<62. (2.7) 

Our choice of log priors is dictated by both physical and statistical considerations. 
From the physical point of view, adopting a prior that is flat in the log of the masses 
means assuming a priori that all orders of magnitudes have the same probability, e.g. that 
the probability of 10 GeV < mo < 100 GeV is the same as the probability 100 GeV < 
TTT-o < 1 TeV. This is appealing because it does not single out any order of magnitude for 



the mass scale in the problem. From the statistical point of view, it was shown in [28| that 
this choice of priors produces a posterior probability that is close to what one would get 
from a naive grid-scan of the likelihood, and maximising the likelihood along the hidden 
dimensions (so-called "profile likelihood"). A detailed statistical analysis shows that the 
average quality of fit of a log prior scan is much better than the average quality of fit 
assuming flat priors ||2^, i.e. priors flat in mo,mi/2- For all those reason, we adopt the 
log prior choice as our default prior for the main results presented in the following. We 
comment on the sensitivity of our results to changes in the prior choice in section |6|. 

2.3 Observables and data 

The experimental values of collider and cosmological observables that we employ to con- 
strain the CMSSM parameter space are listed in Table |^, while the constraints on the 
SM nuisance parameters are given in Table |. We refer to H, |l], H for details about 



the computation of each quantity and for justification of the theoretical errors adopted, 
as well for a detailed description of the likelihood function. In particular, we include the 
measurement of the anomalous magnetic moment of the muon based on e~^e~ data, which 
gives a 3.2a discrepancy with the SM predicted value |43|. As regards BR{B — > Xs^) 



(which has been recently shown to provide an important constraint [^]), for the new SM 
prediction we employ the value of (3.12 it 0.21) x 10~^. We compute the SUSY contribution 



to BR{B Xgj) following the procedure outlined in Refs. |44, 45| which was extended in 



Refs. 47 1 to the case of general flavor mixing. 

Regarding cosmological constraints, we use the determination of the relic abundance 
of cold dark matter based on the 5-year data from WMAP ||4^ to constrain the relic 
abundance ^^^^^ °^ lightest neutralino (assumed to be the sole constituent of dark 
matter). Also, points that do not fulfil the conditions of radiative electroweak symmetry 
breaking and/or give non-physical (tachyonic) solutions are discarded. 



^We also follow the treatment of for the priors on the SM parameters, although this choice does not 
represent an issue as those quantities are directly constrained by experiment and thus do not suffer from 
prior-ambiguity. 
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Observable 



Mean value 



Uncertainties 
a (exper.) r (theor. 



ref. 



Mw 
sin^ 9cs 

BR(B Xs-f) X 10^ 

BR{Bu Tv) X 10^ 



80.398 GeV 
0.23153 
29.5 
3.55 

17.77 ps-i 

1.32 

0.1099 



25 MeV 
16 X 10"^ 
8.8 
0.26 

0.12 ps"i 

0.49 

0.0062 



15 MeV 
15 X 10"^ 
1.0 
0.21 

2.4 ps-i 
0.38 



m 
m 
m 
m 
m 
m 



Limit (95% CL) 



r (theor. 



ref. 



rrih 

rriq 
m-r, 



other sparticle masses 



< 5.8 X 10"^ 

> 114.4 GeV (for an SM-like Higgs) 
/(mO (see ^) 

> 375 GeV 

> 289 GeV 

As in Table 4 of Ref. ptl. 



14% 
3 GeV 
negligible 
5% 
5% 



m 
m 
m 
m 
m 



Table 1: Summary of the observables used in the analysis to constrain the CMSSM parameter 
space. Upper part: Observables for which a positive measurement has been made. (5a^^^^ = 
denotes the discrepancy between the experimental value and the SM prediction of the 



^oxpt 



,SM 



anomalous magnetic moment of the muon [g 
currently exist. 



2)^. Lower part: Observables for which only limits 



SM (nuisance) 
parameter 


Mean value Uncertainty 
/i a (exper.) 


Ref. 


Mt 

l/aem(Mz)*^^ 


172.6 GeV 1.4 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 
127.955 0.03 


|57] 

H] 

|5e] 

m 



Table 2: Experimental mean /i and standard deviation a adopted for the likelihood function for 
SM (nuisance) parameters, assumed to be described by a Gaussian distribution. 



3. Indirect WIMP searches with neutrino telescopes 



As mentioned in the introduction, indirect searches for dark matter are based on the search 
for a neutrino signal from the center of the Earth or the Sun over the known atmospheric 
neutrino background. The most sensitive searches to date are those carried out by neutrino 
telescopes. These are detectors that use huge instrumented volumes of natural ice or water 
both as target and Cherenkov medium to detect the muons^ produced in neutrino-nucleon 
interactions in or around the detector. An array of photomultiplier tubes detects the 
Cherenkov photons with nanosecond resolution and allows for the reconstruction of the 

^Electron or tau-induced particle cascades, in the case of incoming Ve or Vr, do not provide a good enough 
angular resolution for the studies carried out in this paper, so we concentrate here on the i/^ channel only. 
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muon trajectory, and hence of the incoming neutrino direction. The angular resolution 
depends on the incoming neutrino energy, and it is typically a degree at TeV energies. In 
this paper we investigate the signal coming from neutralino annihilations in the center of 
the Sun. 



3.1 Neutrinos from WIMPs annihilation in the Sun 

From an experimental point of view, the relevant observable is the total number of neutrino- 
induced muons from the direction of the Sun per unit time, N^, with an energy above a 
given threshold E^^. This is given by the convolution of the muon effective area of the 
detector, A'^{E^), with the neutrino-induced differential muon flux at its location, 

r+°° (lib 
N, = / dE,Af{E,)-^. (3.1) 

The differential flux of neutrino-induced muons reaching the detector is computed from 
the differential neutrino flux by accounting for both the production cross section of muons 
from neutrinos and their propagation in ice. The differential parent neutrino flux, is 
given by 

dE,~ AttD^ dE,' 

i 

where Fann is the annihilation rate, D is the distance from the source to the detector, the 
sum runs over all WIMP pair annihilation final states and are the branching ratios for 
each channel, each giving rise to the differential neutrino spectrum There is a total of 
29 channels potentially contributing to the sum in Eq. ( |3.2D , although the most important 
ones in the context of the CMSSM (as we show below) involve the final states T^"*"!^", ZZ, 
T^T^ , bb and, when kinematically allowed, tt. The energy range of neutrinos produced by 
WIMP annihilations is bounded by the maximum WIMP mass, which in the case of the 
CMSSM is of order < 1 TeV. 

WIMPs from the galactic halo accrete in the center of the Sun by gravitational capture 
and their number is depleted by annihilation. The capture rate Pcapt on (JsPioc is propor- 
tional to the scattering cross section ag between a WIMP and a nucleus and to the local 
WIMP density, pioc- In order to compute it, one needs to specify both a particle physics 
model (determining (jg) and a halo model, giving the local halo density, pioc, and the ve- 
locity distribution of WIMPs. In this paper we assume an NFW profile with pioc = 0.3 
GeV/cm^ and a Maxwell-Boltzmann velocity distribution with mean v = 270 km/s. As- 
suming a constant capture rate, one can compute the time-dependent density of WIMPs 
in the Sun and from there their annihilation rate, Fanm which is given by 

rann = %^tanh2('i], (3.3) 



where r = l/yT^JaptC^, where Cq oc {(Jav) is proportional to the velocity-averaged annihi- 
lation cross section. For timescales t ^ t, equilibrium between annihilation and capture is 
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established, and thus Tann ~ rcapt/2. This is an important issue for indirect searches since 
only if capture and annihilation are in equilibrium is the neutrino flux from neutralino anni- 
hilations at "full strength". Conversely, if equilibrium is not reached, Fann Tcapt and the 
neutrino signal is correspondingly weaker. In our approach, we compute Fcapt; Ca from our 
theoretical model and galactic halo model, and Tann from Eq. ( |3.3D , assuming t = 4.76 Gyr 
for the age of the solar system. Thus the equilibrium condition can be explicitly checked. 



see section 5.2. 



So far, none of the neutrino telescopes in operation has detected any excess of neutrinos 
from the Sun. This lack of signal has allowed to set limits on the muon flux from neutralino 
annihilations in the Sun, as well as on the annihilation rate of neutralinos. The most 
stringent limits to date are those from the 22-string configuration of IceCube obtained 
with 104 days of livetime during 2007 Q. In the light of this situation, here we focus 
on the expected sensitivity to muons of the final IceCube configuration to describe the 
potential of km^ neutrino telescopes to probe the CMSSM parameter space. 

3.2 The IceCube neutrino telescope 

In this work we will use the proposed original geometry of the IceCube neutrino tele- 



scope 1 49 1 to exemplify the impact that kilometer-scale neutrino detectors can make in 
constraining the CMSSM through searches for an excess in the neutrino flux from the 
Sun's direction. The baseline IceCube geometry consists of 80 strings with 60 digital opti- 
cal modules each, deployed at depths between 1450 m - 2450 m near the geographic South 
Pole. The optical modules are vertically separated by 17 m, while the strings are arranged 
in a triangular grid with an inter-string separation of 125 m. The instrumented volume of 
ice covers 1 km^. Each module contains a 25-cm Hamamatsu photomultiplier tube with 
electronics for in-situ digitization and timing of the photomultiplier signals [^] . 

IceCube is designed to detect the Cherenkov radiation of secondaries produced in 
neutrino interactions in the vicinity of the detector. Although its geometry has been 
optimized to detect ultra-high energy (>TeV) neutrinos from potential cosmic sources 
with sub-degree angular resolution, its energy and angular response are still adequate at 
0(100) GeV to perform searches for dark matter, where the typical neutrino energies are 
of that order or less. A measure of the efficiency of a neutrino telescope to detect muons 
is given by the effective detector area, defined as 



AfiE,) = §^ X (3.4) 



where A'^^ei is the number of events that pass the analysis selection cuts from a test sample 
of A'^gen generated muons that cross a fiducial area A^^'^ . The effective area depends on the 
muon energy, E^, and the incident angle. For the studies carried out in this paper we have 
taken the zenith angle averaged effective area of IceCube as a function of muon energy 
after some simple track quality selection criteria, as published in [49| and shown here in 
Fig. H for completeness (black dots). 
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Figure 1: Effective area of IceCube (final configuration witli 80 strings, from [^) as a function 
of muon energy. Triangles denote trigger level, while black dots (adopted in this work) are the 
effective area after a simple event selection based on track quality parameters. Diamonds and stars 
represent further energy cuts. See |^ for details. 



Recently, an addition of six strings forming a denser core in the middle of the IceCube 
array has been proposed in order to lower the neutrino energy threshold of the detector to 
about 10 GeV [^], the so-called "DeepCore" configuration. This change of the originally 
proposed geometry will significantly increase the sensitivity of IceCube to searches for low 
mass WIMPs. Since performance studies of DeepCore are still under way at the time of 
this writing, we have taken a conservative approach and use the 80-string geometry of 
IceCube as the benchmark for our calculations. 

In order to make statistically sound statements about the capability of IceCube to 
detect or rule out parts of the CMSSM, the expected signal obtained through Eq. ( p?l] ) 
has to be compared with the background consisting, in the ideal detector case, of atmo- 
spheric neutrino-induced events and of an irreducible background of solar atmospheric 
neutrinos. We neglect here any contamination from misreconstructed atmospheric muon 
tracks, assuming the IceCube analysis can reject such tracks during the analysis steps. 
The signal and the background have a different energy and angular dependency, which can 
be exploited in the data analysis to distinguish between them. The signal from WIMP 
annihilations in the Sun is expected to come from a cone of angular size A^' from the 
direction of the Sun, while the atmospheric neutrino background is expected to be, to first 
approximation, isotropic. In this paper we do not employ the signal's energy dependency, 
but we restrict our analysis to a simple counting scenario, assuming that observations of 
atmospheric neutrinos coming from directions away from the Sun can be used to estimate 
the number of background events Nf, above the threshold energy in the signal region. For 
our estimation of the background we follow the results of the simple analysis presented 
in [49 1, according to which ~ 10^ atmospheric neutrinos are expected per year from the 
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whole hemisphere. We wiU here make another conservative choice and take for our our 
analysis a cone of = 3.0° around the Sun position. IceCube is expected to have a bet- 
ter angular resolution than that, especially for tracks of energy 0(100) GeV and above, so 
a dedicated WIMP search analysis with the complete detector is likely to achieve a better 
S/N than what we use here. Rescaled to 1 yr (365 d) of detector live time, we estimate 
about 34 events per year from atmospheric neutrinos from the direction of the Sun. An 
irreducible high-energy neutrino background from the Sun are the neutrinos produced by 
cosmic ray interactions in the corona. Estimations of such effect show that the expected 
number of events in a km^ detector are a few per year [^, ^ . We have used figure 8 in |^ 
to get an estimation of the number of such events that we would expect in IceCube. Since 
the authors assume a flat effective area of 1 km^, independent of energy, we have rescaled 
the prediction one can extract from their figure 8 for an energy threshold of 50 GeV, to a 
more realistic effective area of 0.2 km^ at those energies, based on figure || in this paper. 
This leads to estimating 6 events from the corona per year. Since both the prediction of 
the number of atmospheric neutrinos, as well as the prediction of the number of corona 
events are affected by the uncertainty in the normalization of the cosmic ray spectrum, we 
adopt here a total number of background events from within 3° from the Sun direction in a 
year of data taking of Ni, = 40 yr~^, and we will not specify any further whether they are 
corona or atmospheric events. We neglect in the rest of this paper the uncertainty in the 
estimation of Nf, and we thus assume that its value is perfectly known. We will improve 
on these approximations in a future work. 

We then estimate the signal-to-noise of a WIMP-induced muon rate A'^ over the 
background as 

S/N = -^. (3.5) 



The parameter space region accessible to IceCube is defined by the criterion S/N> 5, i.e., 
as the region where IceCube would obtain an approximate 5cr detection of the WIMP 
signal. This corresponds to an event rate due to WIMPs annihilation > 31.6 yr~^ 
On the other hand, we can also consider the case where a measurement compatible with 
background-only events is obtained. For definiteness, we assume that the observed number 
of events is exactly equal to the predicted background rate, Ni,. We use the Feldman and 
Cousins construction |^2| to derive the 90% upper exclusion limit for the signal event rate, 
NgQ, and we find Nqq = 11.5 yr~^. In other words, in absence of a signal the parameter 
space region predicting > Nqq would be excluded at the 90% CL. When comparing the 
CMSSM predictions with the IceCube reach, we use the S/N> 5 criterium as discovery 
threshold and the > Nqq criterium as the projected exclusion limit. Since the CMSSM 
predicted event rates stretch over more than 10 orders of magnitude (with the majority 
of the parameter space below either the detection or exclusion thresholds), effectively it 
makes almost no difference which one of the two thresholds one adopts to define the reach 
of the detector. For definiteness, when we discuss "the reach" of IceCube we always refer 
to the region delimited by the discovery criterium S/N> 5. 
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4. Indirect detection prospects with IceCube in the CMSSM framework 



We have now assembled all the relevant elements to predict the number of muon events 
in IceCube from WIMPs annihilation in the Sun. Our procedure is as follows. We have 



modified the publicly available SuperBayeS code [33| to link it with DarkSusy 5.0.4 
For every point in the CMSSM parameter space visited by our sampling algorithm, Dark- 
Susy [^] computes the annihilation and capture rates, Fanri! Tcapt, and the branching 
ratios into the different annihilation channels, in Eq. ( |3.2D . These are then used to ob- 
tain the number of neutrino-induced events at the location of the detector by interpolating 
a series of tables with the yields for each annihilation channel obtained via the WimpSim 



code [35, WimpSim employs a Monte Carlo simulation to generate the three-flavour 
neutrino flux at the center of the Sun from WIMP annihilations and to propagate it to 
the Earth including full-flavour oscillations. It has been shown that including the effect of 
neutrino oscillations can have a very substantial impact on the predicted neutrino flux at 
the Earth for some channels [36, 37 [. Here we assume a normal neutrino mass hierarchy 
and the neutrino mixing angle ^13 = 0. Residual uncertainties in the neutrino physics 
(in particular, an inverted mass hierarchy and/or a non-zero ^13 mixing angle |6^) might 
lead to modifications of some of the yields [|6^, although a recent study has not found a 
significant impact [|^. We adopt a muon effective area, A^^, obtained by interpolating the 
values of A'j^{E^) in Fig. Q in the energy range 

£;th = 50 GeV < E <1 TeV, and we use 
it to perform the convolution with the predicted muon spectrum at every point visited in 
the CMSSM parameter space, according to Eq. (3.1). 

Our scanning procedure produces a set of samples {0^^\ ©^^\ . . . } in the CMSSM 
parameter space whose density is proportional to their posterior probability given current 
data, Eq. (^^). For any derived quantity /(O) that is a function of our CMSSM parameters 
(such as the muon flux, or the annihilation and capture rates, computed as described 
above) we can simply plot the distribution of /(O^*-*) {i = 1,...,N) and the latter will 
automatically inherit the statistical properties of the samples, i.e. the distribution of 
derived quantities can also be interpreted in a statistical way. This means that we can 
make probabilistic statements about all sorts of derived quantities (and their correlations) 
that are a function of our CMSSM parameters and/or of any of the derived observables 
(fluxes, rates, etc). 



4.1 Expected number of events and probability of discovery 

We begin by plotting in Fig. |2| the probability distribution of the number of neutrino- 
induced muon events expected for IceCube and as a function of the neutralino mass, after 
all available current constraints given in Tables |l| and ^ have been applied. The posterior 
probability distribution shows at 95% an island with several hundreds of events, above the 
discovery threshold. However, most of the probability density favours region of parameter 
space leading to very few events per year, typically <^ 1. This is highlighted by the right 
panel of Fig. ^, where we plot the ID marginalised distribution for A'^^. The probability 
density peaks at Nfj_ ~ 0.1 yr~^, with a long tail extending until ~ 10"^'' events per year. 
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Figure 2: Left panel: 2D probability distribution for the number of muon events from the the Sun 
for IceCube and neutralino mass in the CMSSM, after all available constraints from current data 
have been applied, with 68% and 95% contours. Right panel: ID probability distribution for the 
number of muon events, with all other parameters marginalised over. The encircled black cross is 
the best fit point. Both plots assume an NFW profile with pioc — 0.3 GeV/cm^ and v = 270 km/s 
for the WIMPs' velocity distribution. 



From this, we find that N^, > 2.3 • lO"'^ at 95% (1-tail limit). Further insight is gained 
by considering the left panel of Fig. |3|, where we plot the expected number of events in 
the {mi/2,mo) plane. We can see that the /i-pole region (thin, vertical sliver at constant 
^1/2) corresponds to the vertical band in Fig. |^ predicting an extremely small number of 
events (and extending all the way down to ~ 10^^'^ yr"^, not shown in the figure). 
This is due to the fact that in the /i-pole region the neutralino mass is low (~ 50 GeV), 
hence the neutrino spectrum is quite soft and this yields muons with energies mostly 
below the IceCube threshold. This exemplifies the importance of lowering the energy 
threshold of neutrino telescopes as much as possible in order to increase the sensitivity 
to all regions of the CMSSM. The stau-coannihilation region (the island at small niQ in 
Fig. ^, corresponding roughly to the 68%, diagonal region in the left panel of Fig. |2|) also 
produces a very small number of events. In this region, the neutralinos are bino-like and 
they annihilate to fermion-anti-fermion pairs, mostly via sfermion exchange. Because the 
neutralino is a Major ana fermion its annihilation cross section to fermion pairs is s-wave 
suppressed and (aav) oc m|/m^. As a consequence, the contribution to the total cross 
section from bb and r+r^ channels is small. Therefore the annihilation rate and capture 
rate are far from being in equilibrium (see section |5.2| for a more thorough discussion), 
and the ensuing neutrino flux is small. A special role is played by the top quark since, if 
this channel is kinematically allowed, one expects a larger cross section and hence a larger 
neutrino flux due to its heavier mass. This requires both a larger value of the neutralino 
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Figure 3: Expected number of neutrino-induced muon events from the Sun for IceCube in the 
('^1/2 1 "lo) plane (left panel, colour coding) and corresponding number of a detection significance 
over the background (right panel). IceCube will be able to probe the parameter space region 
coloured in dark (above the 5cr detection threshold). The density of samples in each plot reflects 
probability density. 



mass and that the neutralino be a mixture state of bino and higgsino. Both conditions 
are satisfied in the focus point region at large {mi/2,mo), which typically corresponds to a 
larger number of events per year, of the order of up to several hundreds. 

We can convert the probability distribution for N^^ into the probability distribution for 
a detection at a given S/N level, estimated according to Eq. ( |3.5D . The distribution of the 
significance level in the (?tT'1/2) '"i-o) plane is shown in the right panel of Fig. |3|, where the 
black points have fluxes above the 5a detection threshold for IceCube. This demonstrate 
that IceCube will be able to probe at high significance the focus point region of the CMSSM. 
We find that IceCube has a 1.7% (3.4%) probability of obtaining a 5a {3a) detection of 
WIMPs from the Sun. However, a large fraction of the CMSSM parameter space favoured 
by current data will remain unaccessible to IceCube. It has to be kept in mind that the 
relative importance of the focus point is difficult to assess with present-day data. In fact, 
a different choice of priors for the CMSSM parameters leads to a significant increase in 
the probability of the focus point, and hence to more optimistic detection prospects for 
IceCube. We comment on the impact of the change of CMSSM priors in section ^, where 
we also compare our results with previous studies. 
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A modification in our astrophysical assumptions and /or in the choice of data employed 
to constrain the CMSSM might lead to modify somewhat the above results. It has been 
suggested ||6^ that the Milky Way might also host a "dark disk" (additionally to the dark 
matter halo) with a population of co-rotating WIMPs whose capture rate for the Sun could 
be up to 1 order of magnitude larger than for halo WIMPs due to their different velocity 
distribution |40|. We can qualitatively estimate the effect of such a dark disk by considering 
an increase in the muon flux by a factor of 10. This would increase the probability of a 
5(7 detection to 4.4%. From the particle physics side, it is worth emphasizing that the 
pull towards the low-mass coannihilation region comes, to a large extent, from one single 
observable, namely the anomalous magnetic moment of the muon, for which a discrepancy 
with the SM value persists at the ~ 3cj level. However, if one removes this somewhat 
controversial observable from the analysis, the relative importance of the focus point region 
(and therefore the detection prospects for IceCube) might increase significantly p^]. 



4.2 Contribution from different annihilation channels 



We now turn to investigating the relative contribution of the different annihilations channels 
to the signal for IceCube. In model-independent analyses of present-day neutrino telescope 
data, it is often assumed that the neutrino spectrum is dominated by one of the 29 available 
final states. Typically, the analysis is carried out for a soft neutrino spectrum, produced 
in quark jets from e.g. bb, and a harder spectrum, ensuing from the decay of r leptons 
or gauge bosons (e.g., ||7^, |^). It is then argued that such choices will roughly cover the 
phenomenology of "typical" points in parameter space. There are two potential issues 
with this procedure. First, once an underlying model is specified (here the CMSSM), the 
relative importance of soft and hard channels will generally vary across the parameter 
space, and in a correlated way. Therefore, limits obtained assuming that only one of 
the channels dominates in the whole parameter space cannot be automatically applied to a 
model-specific scenario. Secondly, even in a purely phenomenological approach, incorrectly 
assuming that one of the channels dominates while the signal really comes from a different 
final state or from a mixture of states leads to potentially large systematic errors in the 
WIMP parameters extracted from indirect detection methods. This has been demonstrated 
in the case of 7 ray observatories and, as we show below, the same caveat applies in 
the context of neutrino telescope data. 



Benchmark points 

Figure ^ shows the differential muon flux for two benchmark points, depicting the contri- 
bution from the most important channels. The left panel (benchmark point A) is for the 
overall best-fit point from our scan, which lies in the coannihilation region, while the right 
panel is for a point in the focus point region, within the reach of IceCube (benchmark point 
B). The parameters values for the two benchmark points are given in Table |3|, along with 
their branching ratios and the fractional contribution to the number of muon events for 
each channel. The branching ratio of benchmark point A is dominated by the bb final state 
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Figure 4: Differential muon flux for the dominant annihilation channels for benchmark point A 
(left panel) and benchmark point B (right panel). The solid, black is the total spectrum, while the 
solid red line depicts the effective area used in the convolution. 

(94.7%), but this channel accounts for only about 8% of the number of events. This is be- 
cause the soft bb spectrum is mostly in the energy region where the effective area of IceCube 
is small, see the left panel of Fig. As a consequence, the harder t^t~ spectrum, whose 
branching ratio is only 1.5%, is responsible for more than half of the expected number of 
events. Overall, this benchmark point will remain outside the reach of IceCube (although 
detection prospects would be more favourable for a DeepCore configuration), as the total 
number of events is = 0.11 per year. Benchmark point B is rather different, as the two 
most important channels in terms of branching ratio are XX ~^ (71%) and final states 
involving gauge bosons (16.6% for W~^W~ and 9.4% for ZZ). After convolution with the 
effective area (see right panel of Fig. Q), the relative importance of those channels changes 
very noticeably, with the final state involving ti contributing about the same fraction of 
events as the final states involving gauge bosons together. 

Systematic error from assuming a single dominating channel 

Let us now consider more in detail benchmark point B, for this represents a fiducial model 
that would be detectable with ~ So" significance by IceCube. If one assumed (as it is done 
in current data analysis from neutrino telescopes) that the spectrum is dominated by one 
single channel (i.e., assuming a branching ratio of 1 for that channel), one would introduce 
a serious systematic error in the analysis. As an illustration, if we assumed that the bb 
was the dominant channel, then we would predict A'^ ~ 0.8 events per year, rather than 
the correct value A^^ = 34.8 yr~^, thus underpredicting the true rate by a factor ~ 40. If 
instead one considered the W'^W^ channel only, then one would obtain A'^^ ~ 64 yr~^, 
thereby overpredicting the number of events by a factor ~ 2. In current practice, the 
observed number of events is used to derive an upper limit to the number of events due 
to a possible WIMP signal. This is then translated into a corresponding limit for the 
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Quantity 


Benchmark A 


Benchmark B 


CMSSM and SM parameters 


TTlQ \\jte V ) 




iZO.D 


zouy.D 


mi/2 (GeV) 




372.9 


666.6 


Ao (GeV) 




728.6 


719.2 


tan P 




18.41 


12.13 


Mt (GeV) 




173.6 


171.0 


mfe(m6)*^ (GeV) 




4.18 


4.34 






0.1179 


0.1170 


l/aem(Mz)^^ 




127.96 


127.90 






Cosmological relic abundance 




0.1145 


0.1022 


Neutralino mass 


(GeV) 


148.2 


267.2 


BR's for dominant channels 


BR(xx ^ hh) 




0.947 


2.1 X 10-3 


BR(xx ^ tt") 







0.717 


BR(xx ^ T+r-) 




0.015 


1.6 X 10-^ 


BR(xx ^ VF+VF- 


-) 


0.011 


0.162 


BR(xx ^ ^^) 




2.5 X 10-3 


0.094 


Muon events for IceCub 


e 


yr-i (total) 




0.11 


34.8 






0.08 


5 X 10-5 









0.46 






0.56 


5 X 10-^ 


Nw+w-/N^ 




0.24 


0.30 






0.06 


0.21 



Table 3: Values of the CMSSM and SM parameters and corresponding branching ratios as well 
as number of muon events in each channel for the two benchmark points discussed in the text. 

Benchmark point A is the overall best-fit for the CMSSM, and it lies in the coannihilation region. 
Benchmark point B is in the focus point region (within ^ 2a of the best-fit) and would be detectable 
with ~ 5cr significance by IceCube. We also give the value of the cosmological relic abundance and 
of the WIMP mass for completeness. 

muon flux, assuming one dominating channel. This step can potentially introduce a large 
systematic error in the flux limit, as the above example illustrates, arising from making 
the extreme assumption of annihilation into only one dominant channel (when the actual 
signal is coming from several channels) . Since the flux constrained by the experiment is the 

product of the branching ratio times the annihilation cross section, this systematic error 
would propagate to the estimation of the annihilation cross section and from there to the 
scattering cross section (derived under the assumption of equilibrium). This is acceptable 
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when data only give an upper limit, as long as one keeps in mind that the limits derived 
under a typical soft (hard) channel are going to be an underestimate (overestimate) of the 
true limits. However, this problem becomes a source of concern once a signal is detected, 
for in that case it is clearly desirable to avoid introducing large systematic errors in the 
estimates of the annihilation cross section and of the neutralino mass (the latter effect 
coming from employing the wrong spectral shape when one assumes a single dominating 
channel). A methodology that accounts for this and avoids introducing systematic biases 
will be presented elsewhere. 

In general, we can introduce a factor f^^^^ which measures the error introduced in the 
annihilation and scattering cross section estimation when assuming that a single channel i 
dominates the annihilation (assigning it a branching ratio of unity) 

where Nj represents the number of events due to annihilation channel i and BR(xx ~^ ^) 
is the true branching ratio for that final state. The systematic error is largest for points 
in parameter space whose contribution to the fractional number of events is very different 
from their branching ratio. Notice that f^^^^ ^ 1 when the branching ratio is close to unity 
but Ni/N^ <C 1, which is the case for the soft channels, while ff^^^ <C 1 when the signal 
is dominated by one or more channels with a small branching ratio (typically the case for 
i = W+,W-,ZZ OTT+T-). Fig. I shows the value of /f* for i = hb,T+T~,W+W- and « 
in the plane of the neutralino mass and the muon flux^. The figure only shows samples in 
the CMSSM parameter space in the IceCube discovery region. The final state which leads 
to the least bias is XX ~^ W~^W~ (less than 50% error across most of the accessible region), 
while assuming a bb {tt) final state leads to a systematic underestimation (overestimation) 
of order ~ 100 or larger. 



Final-state contributions to the signal 

We now turn to investigating the relative importance of the various final states across the 
whole CMSSM parameter space favoured by current data. The most important annihilation 
channels are summarized in Table |^, while in Fig. ^ we plot the branching ratios (top row) 
and the fractional number of muon events (bottom row) in the ("ii/2)^o) plane for 4 out 
of the 5 annihilation channels considered in the table (the plots for ZZ are qualitatively 
similar to the ones involving W~^W~ and are not shown). The same quantities but in 
the (m^,<I>^) plane are shown in Fig. ^. Although the dominant annihilation mode is 
XX ~^ bb, whose average branching ratio is 0.67, this soft channel actually only contributes 
an average 4% to the total number of events. As argued above, this is because most of the 
differential muon flux from bb is at low energies ( ^ 10 GeV, see Fig. below the threshold 
energy adopted for IceCube. Therefore after convolution with A'^^ the resulting number 
of muon events is a small fraction of the total. Clearly, a lowering of the threshold energy 

''For the 66 final state there are several points where this channel contributes no events at all. Such 
points have been arbitrarily mapped to a value /ff''* = 100. 
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Figure 5: Systematic error factor induced in the reconstructed annihilation rate from (incorrectly) 
assuming that the signal is dominated exclusively by one of the four channels shown (from left to 
right: bb , r+T~, W^W~ and ti)^ as a function of the muon flux above 1 GeV and the neutralino 
mass (notice that the muon flux has been converted to the conventional 1 GcV threshold, and that 
the scale of the vertical axis is linear). Green points are the least biased, having a systematic error 
below a factor of 3. We show only samples within the discovery region of IceCube. 



Channel i 


(BR)cMSSM {Ni/N^)cMS^M 


(BR)lccCube (-/Vi/-/V;i) IceCube 


XX 

XX T+T- 
XX it 

XX ^ w+w- 
XX zz 


0.67 0.04 
0.13 0.65 
0.07 0.08 
0.02 0.11 
5 X 10-3 0.03 


6.6 X 10-2 6.7 X 10-'' 
8.6 X 10-3 3.0 X 10-2 
0.26 0.14 
0.49 0.60 
0.14 0.20 



Table 4: Dominant annihilation channels and the associated fractional number of events, aver- 
aged over the posterior of the CMSSM parameter space (columns 2 and 3) and averaged over the 
parameter space accessible to IceCube (defined as the region with S/N > 5, columns 4 and 5). 



and an increase of the effective area at small muon energies would allow the DeepCore 
configuration to increase significantly its sensitivity to this soft channel. The second most 
important channel (in terms of its average branching ratio) is XX ~^ t~^t~ , which leads to a 
harder spectrum and therefore to a larger fraction of the number of events after convolution 
with the effective area. Indeed, the t~^t~ is responsible for ~ 65% of the total number 
of events. By comparing the plot of the fractional number of events in the t^t^ channel 
(second panel from the left in the bottom row of Fig. |^ with the right panel of Fig. |3|, 
it is clear that this annihilation mode dominates in a region of the CMSSM parameter 
space which will remain undetectable by IceCube. On the other hand, if we consider the 
parameter space region within the reach of IceCube, we see from Fig. ^ that the most 
relevant channels will be the ones with gauge bosons final states (VF~^VF- and, not shown, 
ZZ), as well as the tt final state, for the reasons explained above. In the parameter space 
accessible to IceCube, the two final states involving gauge bosons account together for 
about 80% of the signal on average, while the ti channel accounts for an average 14%, see 
Table § 

In summary, in the context of the CMSSM the overall dominating annihilation mode 
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Figure 6: Top row: branching ratios in the ('Ti]^/27"^o) plane for the most important annihilation 
channels in the CMSSM, from left to right: bh, t^t^ , W^W~ and tt (the plots for ZZ are similar to 
the ones for W^W~ and are not shown). Bottom row: fractional number of muon events expected 
in IceCube from each channel. 



in terms of its branching ratio is XX ~^ bb, although we have found that this contributes 
only a very small fraction to the number of events that are within reach of IceCube. In 
the parameter space region accessible to IceCube, though, the most important annihilation 
channels (which dominate the number of events) are the ones involving gauge bosons final 
states as well as xX ~^ it- We have also found that the least biased choice for a single- 
channel reconstruction is xX ~^ W~^W~ , which together with xX ~^ (which has a 
very similar spectral shape) accounts for ~ 80% of the number of events in the CMSSM 
parameter space accessible to IceCube. 



5. Comparison with direct detection techniques 



5.1 Ruling out the CMSSM with IceCube and direct detection 

It is interesting to consider the cross-correlation between the signal in a neutrino telescope 
and the expected signal from direct detection methods, as this gives a consistency check 
of the model. To this aim, we computed the neutralino spin-independent scattering cross 
section off a proton according to the prescription given in Ref. |71|. When comparing the 
indirect detection prospects using neutrino telescopes with direct detection experiments, 
it is important to keep in mind that predictions for the latter also suffer from a number of 
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Figure 7: As in Fig. ^but in the (m^, <i>^) plane instead, where is the muon flux above 1 GeV. 



sizable uncertainties. The value of the neutralino-nucleon elastic scattering cross section 
is subject to up to a factor of ~ 3 (~ 10) uncertainty at 68% C.L. (at 95% C.L.) due to 
the estimate of the hadronic matrix elements |72|. On the other hand, uncertainties in 
the form factors appear to be subdominant |73]. When comparing experimental counts 
with theoretical predictions, the local WIMP distribution has to be taken into account, 
and in particular the local WIMP density pioc and their velocity distribution. Although 
it has generically been assumed that pioc suffers from uncertainties of up to a factor of 
2 or 3, recent studies have demonstrated that this crucial quantity can be determined 
to 10%-20% accuracy by using a variety of kinematic data to constrain a parameterized 
model of the galactic halo |74]. Furthermore, given at least two experiments with different 
target nuclei masses, the WIMP mass and scattering cross section can be reconstructed 
independently of the velocity distribution of the halo [75|. It is presently difficult to quantify 
the potential impact of other astrophysical sources of uncertainties, such as for example a 
non-smooth component to the local density and the effect of a co-rotating dark matter halo. 
For definiteness, when showing below direct detection constraints on the spin-independent 
cross section from direct detection experiments, we adhere to the standard convention of 
assuming a smooth halo with p\oc = 0.3 GeV/cm'^ and a Maxwellian velocity distribution 
with mean velocity v = 220 km/s. 

Figure § shows the 68% and 95% favoured regions in the plane spanned by the number 
of muon events expected for IceCube vs the neutralino scattering cross section off a nu- 
cleon (spin-independent scattering off a proton, a^, in the left panel and spin-dependent 
scattering off a neutron, <t^^, in the right panel). As expected, the region of parame- 
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ter space leading to a large number of events for IceCube also predicts larger scattering 
cross sections. Current 90% upper limits for the spin-independent cross section are around 
< 5 X 10-8 pb ^ ^, |7| for a 100 GeV neutralino, while the spin-dependent cross sec- 
tions are constrained by direct detection experiment ^ 2 x 10^^ pb and af^^ ^ 10^^ pb, 
for m-^ ^ 100 GeV 81, 82, 83|. The SuperKamiokande neutrino detector puts more 



stringent constraints on the spin dependent scattering cross section off a proton, with an 
upper limit C7p^ ^ 4 x 10~^ pb |84], while the current IceCube data provide an upper limit 



(Tp^ ^ 3 X 10 ^ pb (although this limit applies for a slightly higher neutralino mass of 250 
GeV) i. 

It follows from the left panel of Fig. ^ that should IceCube detect neutrinos from 
the Sun at high significance, then the corresponding spin-independent cross section must 
satisfy 0"^^ ^ 0.9 x 10^^ pb. Therefore, another factor of ~ 5 improvement in the direct 
detection constraints on the spin-independent scattering cross section would rule out the 
CMSSM parameter space region in which IceCube would be able to detect a signal. On 
the other hand, if IceCube does detect neutrinos from the Sun at high significance while 
direct detection experiments fail to find a signal above cXp^ > 0.9 x 10"^ pb, this would 
strongly disfavour the CMSSM, at least given the conventional astrophysical assumptions 
made in the present analysis. Changes in the local WIMP density would equally affect the 
signal for both direct detection and neutrino telescopes, at least for points satisfying the 
equilibrium condition, so the relative power would remain unchanged (although of course 
the absolute prospects of detection would be modified). The same applies with respect 
to the uncertainties in computing the neutralino-nucleon scattering cross section. In the 
presence of a dark disk, the parameter space region accessible to IceCube increases as 
discussed above, but it still corresponds to values of the spin-independent scattering cross 
section > 0.7 x 10"^ pb. We have found that the above results are robust within 20% 
under a change of priors for the CMSSM parameters, as discussed in section 

5.2 The capture annihilation equilibrium condition 

In comparing the constraints obtained from searches using neutrino telescopes with the 
limits obtained via direct detection techniques, it is often assumed that equilibrium is 
reached between the annihilation and capture rates in the Sun, (see e.g., and for a 
recent discussion [^). This means from Eq. (|3.3| ) that Fann = rcapt/2, and therefore lim- 
its on Tann from neutrino telescopes can be directly translated into corresponding limits 
for the scattering cross sections (either spin-dependent or spin-independent), which are 
proportional to Fcapt- It is thus important to verify whether the equilibrium condition is 
indeed achieved across the CMSSM parameter space. We plot the 2D probability distribu- 
tion for Fann and Fcapt in Fig. ^ The region of the CMSSM parameter space that indeed 
has reached equilibrium corresponds to the focus point (around the dashed, diagonal green 
line, in the region with larger values of Fann); and to parts of the stau coannihilation region 



(cf. Fig. 10, showing the equilibrium condition in the (mi/2;^o) plane). However, we can 
see that in other parts of the coannihilation region (branching off diagonally from the best 
fit point, denoted by an encircled cross) and the /i-pole region (isolated diagonal island 
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Figure 8: 2D probability distribution for the number of muon events in IceCube vs the spin- 
independent (left panel) and spin-dependent (right panel) scattering cross section (the theoretical 
predictions are very similar for scattering off a proton or off a neutron). The vertical lines give 
the current best 90% upper limits from direct detection experiments ||7^, |7^ (for = 100 GeV), 
the horizontal line delimits the parameter space that IceCube will be able to exclude with 90% 
confidence. 



at 95% probability) the annihilation rate is very significantly lower than the capture rate, 
and hence equilibrium is not attained. This also explains the low number of muon events 
expected in those parts of parameter space. 

In order to quantify the extent to which the CMSSM parameter space is in equilibrium, 
we define the equilibrium parameter 

^ = log^, (5.1) 

J- capt 

where S' = denotes equilibrium. We consider that equilibrium is not reached if Fann 
is below a factor of 2 the equilibrium condition, i.e. for < —0.3. The value of the 



equilibrium parameter in the {1^1/2, mo) plane is plotted in Fig. 10, confirming that while 
points in the focus point have reached equilibrium, most samples in the coannhiliation 
region and in the h-pole region are out of equilibrium. Indeed, ~ 66% of the samples across 
the parameter space violate the equilibrium condition. However, the region accessible to 
IceCube is delimited by Fcapt ^ 10^^ s~^ and Fann ^ 5x 10^^ s~^, and this part of parameter 
space is indeed in equilibrium, as would be expected. 

6. Prior dependency 

The dominant source of uncertainty in our prediction is the choice of priors for the CMSSM 
scalar and gaugino masses. All results presented until here employed a log prior, but other 
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Figure 9: 2D probability distribution for the annihilation rate in the center of the Sun, Fann, vs 
the capture rate of WIMPs in the Sun, Fcapt, in the CMSSM. The diagonal, dashed green line 
denotes the equilibrium condition, where capture and annihilation are balanced. 



choices are possible (for more details, see e.g. [27, 28, Q). Since the prior choice sets the 
probabilistic measure of for the parameter space, it has been shown that certain classes of 
priors can naturally quantify the physical fine-tuning of the theory |35, 36]. One however 
expects the choice of prior to become irrelevant once the data included in the likelihood 
are strongly constraining, a limit not yet attained by present-day data but which will be 
obtained by future LHC measurements |42|. The main effect of a change of priors is to 
modify the relative importance of the focus point region with respect to the coannihilation 
region. Rather than performing a complete exploration of several classes of priors, here we 
simply investigate the robustness of our results with respect to a change of prior on the 
CMSSM masses by showing in Fig. 11 results obtained assuming a flat prior on (mo, ^1/2) 
over the range 50 GeV < 1710,1711/2 < 4 TeV, with all other aspects of the scan unchanged 
(compare with Fig. adopting a log prior). Such a choice is expected to roughly bracket 
the results one would obtain when analysis a larger class of priors. We can see that under 
a flat prior the preference for the low mass coannhilation region is attenuated, and there 
is a much larger statistical weight on the focus point region (which will be accessible to 
IceCube). This comes about because under the flat prior, the volume encompassed by the 
mass range between e.g. 1 TeV and 4 TeV is a factor 100 larger than between 100 GeV 
and 400 GeV, while under the log prior the ratio of the volumes of the two regions is unity. 
Since the preference for the low-mass region encoded in the likelihood is not sufficiently 
strong to override the much larger prior volume of the focus point under the fiat prior, the 
ensuing constraints remain at this point somewhat prior dependent. 

Using the fiat prior, we find that the parameter space region accessible to IceCube at 
the 5(T level is larger than what we had found under the log prior, and is now about 12.5% 
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Figure 10: The distribution of the equihbrium parameter (a defined in Eq. (5.1) (colour coding) 
in the (mi/2,'7io) plane. Regions in parameter space with (a > —0.3 (dark points) have reached 
equilibrium between capture and annihilation in the center of the Sun. Although about 66% of 
the parameter violates the equilibrium condition, the focus point region accessible to IceCube is in 
equilibrium. 



of the total parameter space. An even more pronounced boost of the detection prospects is 
obtained if one introduces theoretical priors that further favour the focus point region, as 
has been done in Ref. [^]. A detailed comparison with the rather more optimistic detection 
prospects for IceCube claimed in Ref. |27] is difficult, as there can be several differences 
in our analyses. In Ref. [^] not many details are given about their treatment of neutrino 
oscillations, nor about the effective area employed for the telescope. In our case, we have 
used a state-of-the-art code for our prediction of the neutrino flux at the Earth, which 
fully accounts for neutrino oscillations, and convolved the flux with the energy-dependent 
effective area obtained by Monte Carlo studies of the telescope. Although their "flat tan (3 
prior" corresponds to our flat prior analysis, they seem to obtain a sizable probability of 
several thousands of events per year (compare the dashed line in the right panel of Fig. 8 



in Ref. with the right panel of our Fig. 11). 

The authors in Ref. have studied the effect of a possible dark matter disk in the 
galaxy on the capture of neutralinos in the Sun and the corresponding expected muon flux 
at the Earth. They, however, do not make concrete predictions on the number of expected 
events in any particular detector, as we do in this work, but just compare the predicted flux 
with existing experimental limits and with projected IceCube limits, both derived using 
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Figure 11: Left panel: probability distribution for the number of muon events vs neutralino 
mass, adopting a flat prior on the CMSSM parameters. Right panel: marginalized ID probability 
distribution for the number of muon events. The black encircled cross is the best fit (cf. Fig. ^, for 
a log prior). 



the assumptions that a hard channel dominates the spectrum (which we have shown can 
introduce systematic errors). If we, very simplistically, mimic the effect of such disk by 
increasing our predictions by an order of magnitude (the effect found in Figure 1 in pC|]) 



we obtain very similar results to Ref. |40| in terms of the flux expected in IceCube. Using 
flat priors and in the presence of a dark disk, the IceCube accessible region increases to 
about 42% of the CMSSM parameter space. This boost in the probability of a discovery 
comes about because under the flat prior there is a large region of high probability density 
coming from the focus point that sits just below the IceCube sensitivity under normal 



astrophysical assumptions (see the left panel of Fig. 11), but when a dark disk is added 
to the picture, its larger capture rate decisively thrusts this region above the detection 
threshold. However, this scenario has to be considered as very optimistic, in that it relies 
both on astrophysical and statistical effects (coming from theoretical prior consideration) 
to increase by a large factor the probability of discovery. Furthermore, it would seem from 
the left panel of Figure 4 in Ref. |^] that such a large boost to the muon flux would be 
already under pressure from current neutrino telescope data, and therefore if the dark disk 
exists its properties remain at the moment quite uncertain. In conclusion, it seems more 
reasonable to consider the ~ 5% probability of discovery obtained under the log prior and 
the dark disk assumption as a more realistic upper limit on the IceCube discovery potential. 

7. Conclusions 

We have studied in detail the ability of the nominal IceCube configuration to probe the 
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Figure 12: Signal-to-noise levels over the background for the WIMP signal from the Sun (cf. with 
Fig. H, which adopts a log prior on the CMSSM masses instead). 



favoured parameter space of the CMSSM, including a realistic description of the detector 
and the expected background from atmospheric neutrinos and the background from the 
solar corona. We have however taken several conservative assumptions on the performance 
of IceCube in our calculations, so the results presented in this paper are conservative. The 
detection prospects of a WIMP induced signal from the Sun are best for IceCube in the 
focus point region, where equilibrium between capture and annihilation rate in the Sun is 
achieved and the signal is correspondingly boosted. We notice that this important part of 
parameter space will not be accessible to the LHC. Depending on the details of the analysis 
and on the astrophysical assumptions, we found that IceCube has between 2% and 12% 
probability of detecting a WIMP signal at more than 5cr significance. 

We have shown that the assumption that the signal is dominated by only one annihila- 
tion channel, as has become common practice in several experimental analyses, can lead to 
a considerable systematic bias in the conversion to the WIMP scattering cross section. The 
soft channels such as XX ~^ ^'I'e the most prone to this effect, while final states involving 
gauge bosons are the least affected. We therefore recommend adopting the xX ~^ W~^W~ 
channel for the least bias single-channel reconstruction in the context of the CMSSM. 

By comparing the expected signal for IceCube with direct detection prospects, we have 
concluded that the CMSSM will be ruled out if IceCube does detect a WIMP signal from the 
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Sun while direct detection experiments do not find anything above ^ 0.7 x 10"*^ pb (for 
the more conservative choice of astrophysical model). We further find that increasing the 
sensitivity of IceCube to lower muon energies will boost its discovery potential of neutralino 
dark matter, since it will improve access to the XX ~^ annihilation channel, which is 
the dominating mode in the CMSSM. However this mode produces mainly low energy 
neutrinos, below the threshold of the original IceCube configuration. In this respect, the 
proposed low-energy DeepCore extension of IceCube will be an ideal instrument to focus 
on the relevant CMSSM areas of neutralino parameter space. 
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